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In designing linac structures for multi-bunch applications we are often in- 
terested in estimating the effect of relatively weak multi-bunch beam break- 
up (BBU), due to the somewhat complicated wakefields of detuned struc- 
tures. This, for example, is the case for the injector linacs of the JLC/NLC 
linear collider project (see Ref. |l|]). Deriving an analytical formula for such 
a problem is the subject of this report. Note that the more studied multi- 
bunch BBU problem, i.e. the effect on a bunch train of a single strong mode, 
the so-called "cumulative beam break-up instability" (see, e.g. Ref. 0), is a 
somewhat different problem, and one for which the approach presented here 
is probably not very useful. 

In Ref. [|| an analytical formula for single-bunch beam break-up in a 
smooth focusing linac, for the case without energy spread in the beam, is 
derived, the so-called Chao-Richter-Yao (CRY) model for beam break-up. 
Suppose the beam is initially offset from the accelerator axis. The beam 
break-up downstream is characterized by a strength parameter T(i, s), where 
t represents position within the bunch, and s position along the linac. When 
T(t, s) is small compared to 1, the growth in betatron amplitude in the linac 
is proportional to this parameter. When applied to the special case of a 
uniform longitudinal charge distribution, and a linearly growing wakefield, 
the result of the calculation becomes especially simple. In this case the 
growth in orbit amplitude is given as an asymptotic power series in T(t, s), 
and the series can be summed to give a closed form, asymptotic solution for 
single-bunch BBU. The derivation of an analytic formula for multi-bunch 
BBU is almost a trivial modification of the CRY formalism. We will here 
reproduce the important features of the single-bunch derivation of Ref. |J 
(with slightly modified notation), and then show how it can be modified to 
obtain a result applicable to multi-bunch BBU. 

Let us consider the case of single-bunch beam break-up, where a beam 
is initially offset by distance yo in a linac with acceleration and smooth 
focusing. We assume that there is no energy spread within the beam. The 
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equation of motion is 
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with y(t, s) the bunch offset, a function of position within the bunch t, and 
position along the linac s; with E the beam energy, [l//3(s)] the betatron 
wave number, eNt the total bunch charge, \t{t) the longitudinal charge 
distribution, and W(t) the short-range dipole wakefield. Our convention is 
that negative values of t are toward the front of the bunch. Let us, for the 
moment, limit ourselves to the problem of no acceleration and /3 a constant. 
A. Chao in Ref. Q expands the solution to the equation of motion for this 
problem in a perturbation series 

oo 
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with the first term given by free betatron oscillation [y^ = yo cos(s//3)]. He 
then shows that the solution for the higher terms at position s = L, after 
many betatron oscillations, is given by 
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with 
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and Rq(z) = 1. An observable y is meant to be the real part of Eq. ||. 
The effects of adiabatic acceleration, i.e. sufficiently slow acceleration so 
that the energy doubling distance is large compared to the betatron wave 
length, and f3 not constant, can be added by simply replacing (f3/E) in 
Eq. H by (j3/E), where angle brackets indicate averaging along the linac 
from s = to s = L.Q For example, if the lattice is such that (3 ~ E^ 

1 Note that the terms yoe lL ^ 13 in Eq. ||, related to free betatron oscillation, also need 
to be modified in well-known ways to reflect the dependence of /3 on E. It is the other 
terms, however, which characterize BBU, in which we are interested. 
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then (/3/E) = ((3 /E )g{E f /E , (), where subscripts "0" and "/" signify, 
respectively, initial and final parameters, and 

9M) = 1(4 Z t) \P~Et\. (5) 
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Chao then shows that for certain simple combinations of bunch shape and 
wake function shape the integrals in Eq. |] can be performed analytically, and 
the result becomes an asymptotic series in powers of a strength parameter. 
For example, for the case of a uniform charge distribution of length i (with 
the front of the bunch at t = 0), and a wake that varies as W = W't, the 
strength parameter is 

e 2 N t LW't 2 f3 

Mt,L) = g{E f /E ,Q ■ (6) 

If T is small compared to 1, the growth is well approximated by T. If T 
is large, the sum over all terms can be performed to give a closed form, 
asymptotic expression. 

For multi-bunch BBU we are mainly concerned with the interaction of 
the different bunches in the train, and will ignore wakefield forces within 
bunches. The derivation is nearly identical to that for the single-bunch 
BBU. However, in the equation of motion, Eq. |], the independent variable 
t is no longer a continuous variable, but rather t takes on discrete values 
t m = mAt, where m is a bunch index and At is the bunch spacing. Also, W 
now represents the long-range wakefield. Let us assume that there are M, 
equally populated bunches in a train; i.e. Nt = MN, with N the particles 
per bunch. The solution is again expanded in a perturbation series. In the 
solution, Eq. ||, the R n (t), which are smooth functions of t, are replaced by 

_. m— 1 

(with 72.^ = 1), which is a function of a discrete parameter, the bunch index 

m. Note that TZ^ = S m /M, with S m the sum wake. 

Generally the sums in Eq. |7] cannot be given in closed form, and there- 
fore a closed, asymptotic expression for multi-bunch BBU cannot be given. 
We can still, however, numerically compute the individual terms equivalent 
to Eq. | for the single bunch case. For example, the first order term in 
amplitude growth is given by 

Tm= e2N ^ mP ° 9(Ef/Eo,0 [m = l,...,M]. (8) 
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If this term is small compared to 1 for all m, then BBU is well characterized 
by T. If it is not small, though not extremely large, the next higher terms can 
be computed and their contribution added. For T very large, this approach 
may not be very useful. 

From our derivation we see that there is nothing that fundamentally 
distinguishes our BBU solution from a single-bunch BBU solution. If we 
consider again the single-bunch calculation, for the case of a uniform charge 
distribution of length £, we see that we need to perform the integrations 
for R n in Eq. ||. If we do the integrations numerically, by dividing the 
integrals into discrete steps t n = (n — l)At and then performing quadrature 
by rectangular rule, we end up with Eq. [?] with M = t/At. The solution is 
the same as our multi-bunch solution. What distinguishes the multi-bunch 
from the single-bunch problem is that the wakefield for the multi-bunch case 
is not normally monotonic and does not vary smoothly with longitudinal 
position. For such it may be more difficult to decide how many terms 

are needed for the sum to converge. 

In Fig. H we give a numerical example: the NLC prelinac with the op- 
timized 3-7r/4 S-band structure, but with 10 -5 systematic frequency errors, 
with the nominal (2.8 ns) bunch spacing (see Ref. The diamonds give 
the first order (a) and the second order (b) perturbation terms. The crosses 
in (a) give the results of a smooth focusing simulation program (taking 
f3 ~ E 1 / 2 ), where the free betatron term has been removed. We see that the 
agreement is very good; i.e. the first order term is a good approximation 
to the simulation results. In (b) we note that the next order term is much 
smaller. 
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Figure 1: A numerical example: the NLC prelinac with the optimized 
37r/4 S-band structure, but with 10~ 5 systematic frequency errors, with 
the nominal (2.8 ns) bunch spacing (see Ref. Q). Here iV = 1.2 x 10 10 , 
E Q = 1.98 GeV, Ef = 10. GeV, L = 558 m; the rms of the sum wake 
Srms = -005 MV/nC/m 2 . The diamonds give the first order (a) and the 
second order (b) perturbation terms. The crosses in (a) give smooth focusing 
simulation results with the free betatron term removed. 
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